Matlab code for Model pre-A

function p_dprime= model_pre_a(t,p_d)
%% Units
% [p_dprime] = [mol*dm^-3*s^-1]
% [s_d] = [mol*dm^-3*s^-1]
% [d_d] = [s^-1]
% [p_d] = [mol*dm^-3]

%% Constants
s_d = 2.4998*1e-0;      % production rate of dioxygenase by transcription and translation
d_d = 0.000289;         % degradation rate of dioxygenase


%% Equation
p_dprime = s_d - d_d*p_d;     

end


t=0:1:5000;
[t,p] = ode45(@model_pre_a,t,0);

figure
plot(t,p)